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the surface of the Temsheet were used to confirm the presence 
of the Taylor-Gort ler vortices . Average Nusselt number 
measurements were used to investigate their effect on heat 
transfer from the curved test section. Experiments were 
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mically developed and thermally developing. It was confirmed 
that the Taylor-Gortler vortices enhanced the heat transfer 
process in the curved section. The experimental data was 
compared with several analytic studies with good results 
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ABSTRACT 



Experimental investigations have been conducted to 
examine the general problem of convective heat transfer in a 
rectangular cross section channel with both straight and 
curved passages. The experimental configuration was modeled 
as infinite parallel plates with one wall at a constant heat 
flux and the opposite wall adiabatic. The air flow through 
the channel was heated through the Joulean heating of a wall 
composed of Temsheet. Liquid crystals applied to the surface 
of the Temsheet were used to confirm the presence of the 
Taylor-Gortler vortices. Average Nusselt number measurements 
were used to investigate their effect on heat transfer from 
the curved test section. Experiments were conducted for 
steady state laminar flow that was hydrodynamically developed 
and thermally developing . It was confirmed that the Taylor- 
Gortler vortices enhanced the heat transfer process in the 
curved section. The experimental data was compared with 
several analytic studies with good results . 
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Dean number 
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atmospheric pressure N/m 
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Meaning Units 

heat transferred by radiation W 

heat transferred from the inner wall W 

heat transferred from the outer wall W 

gas constant for air J/Kg °K 
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resistor 
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total radiation resistance m 
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bulk temperature of the flow °C 
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average flow outlet temperature after °C 

each test section 

average temperature of inner wall °K 

average temperature of outer wall °K 

voltage across precision resistor V 

voltage across wall heater (Temsheet ) V 

dimensionless axial length coordinate 
emissivity of inner wall (Plexiglas) 
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I. INTRODUCTION 



A. DESCRIPTION OF TAYLOR-GORTLER VORTICES 

A considerable amount of data and indirect evidence, [l, 

2, 3] among others, has shown that the fully developed laminar 
flow of viscous fluids along a concave wall does not remain 
two dimensional. Instead, the flow forms spiral vortices in 
counter-rotating pairs which frequently possess a regularly 
spaced cellular structure. This phenomenon, known as 
Taylor-Gortler vortices, occurs as a result of flow insta- 
bility induced by variations in the centrifugal forces acting 
on fluid particles. The Taylor-Gortler vortices are laminar 
in nature with their axes aligned in the direction of the 
principal fluid flow and with secondary flow velocities in 
both perpendicular directions . An example of this flow 
pattern is shown in Figure 1. 

In a channel that is curved in the streamwise direction, 
such as that used in this study, the fluid particles near 
the center are subjected to higher centrifugal forces than 
those experienced by the slower moving fluid particles near 
the boundary of the concave wall. Consequently, the tendency 
is for the fluid particles in the center of the channel to 
move outward to the concave wall. Under this influence, 
particles near the wall are displaced initially in the 
spanwise direction and subsequently continue radially inward, 
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Figure 1. Taylor-Gortler Vortex Flow Pattern 
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replacing the central fluid particles. Once in the center 
of the channel the fluid particle will come under the in- 
fluence of the higher free stream velocity and resultant 
centrifugal forces. This action forces movement back toward 
the concave wall causing the entire cyclic rotation process 
to propagate . 

Although the Taylor-Gortler vortices are laminar, it is 
quite possible that they play an important role in the transi- 
tion from the laminar to the turbulent flow regime. [2] 
Additionally, the presence of the vortices have been thought 
to account for the increase in heat transfer from curved 
walls due to mass transfer. [4] 

The Taylor-Gortler vortices also have similarities to 
other vortex flow patterns such as the longitudinal vortex 
rolls developed in the laminar forced convection heating 
of fluid layers between parallel plates. [5] In another 
field it has further been hypothesized that the cross-hatching 
phenomenon observed in re-entry vehicles can be explained 
by the presence of centrifugally induced secondary motion 
characteristic of the Taylor-Gortler vortices. [6] 

It is apparent that the understanding of these vortices 
and their effect on heat transfer and fluid flow character- 
istics could significantly effect many areas of fluid 
mechanics. Recently it has been proposed [7, 8, 9j that 
concave wall curvature can have an advantageous effect on 
film cooling over turbine blades . 
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Many applications of fluid mechanics and heat transfer 
involve curved surfaces or curved channels and the under- 
standing of the fluid mechanics in this geometry will greatly 
enhance engineering designs and capabilities. 

B. BRIEF HISTORY 

Instability of an inviscid fluid in a curved flow was 
first considered by Lord Raleigh in 1916. [10 ] He concluded 

that the stability criterion for an inviscid fluid had as a 
necessary condition that the circulation always increased 
with the radius. In essence, this meant that the product 
of the local circumferential velocity and the corresponding 
local radius of curvature cannot decrease with increasing 
radial distance. 

In 1923, G. I. Taylor expanded the scope of Raleigh’s 
stability criterion to include viscous fluids. Cl, 11] His 
extensive analytic and experimental investigation focused on 
the flow of viscous fluids between a rotating inner cylin- 
der and a stationary outer cylinder. Taylor ascertained that 
such couette flows become unstable when a characteristic, 
non-dimensional parameter exceeded a critical value of 41. 
Referred to as the Taylor number, this parameter was defined 
as : 



Ta 




i 



l 
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In this relationship d represents the spacing between the 
concentric cylinders, assumed small when compared to the 
radius of curvature, R^ . Re is the Reynolds number based on 
U, the circumferential velocity of the inner cylinder and d, 
the annulus width. When flow conditions were such that the 
Taylor number exceeded this critical value, the secondary 
flow observed in the gap between the stationary outer cylin- 
der and the concentrically rotating inner cylinder was in the 
form of longitudinal vortices whose axes were aligned in the 
circumferential direction. These vortices were given the name 
Taylor vortices. 

A similar form of instability occurs when a viscous 
fluid flows in a curved channel due to a pressure gradient 
acting along the concave wall. W. R. Dean [ 12 ] first addressed 
this problem analytically in 1923. Specifically, he consi- 
dered flow in a channel formed by two concentric cylinders with 
the restriction that the cylinder gap distance be small com- 
pared to the radius of the inner cylinder. Dean concluded 
that the flow instability would initiate, and the secondary 
flow vortices (similar to the Taylor vortices) would form, 
when the Dean number was greater than 36. The Dean number 
was then defined as: 

De = Re / 
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In this case, the Reynolds number was based on the mean velo- 
city of the undisturbed flow and d represented the channel 
half width. The analytical work by Dean was later verified 
by W. H. Reid [13] using an approximate numerical solution. 

In 1940, H. Gortler [2] studied the influence of small 
disturbances on the stability of laminar boundary layer 
profiles. Gortler found that these disturbances closely 
resembled those described by G . I. Taylor which led to fluid 
flow instability in the form of vortices. In his numerical 
calculations, Gortler concluded that the amplified distur- 
bances were produced only on the concave walls; nonetheless, 
the overall flow regime seemed to remain laminar. The appro- 
ximate analytic results that Gortler had developed were 
initially verified with an exact solution by G. Hammerlin as 
reported by H. Schlicting [14], and later substantiated by 
an extensive numerical analysis by A. M. 0. Smith [3]. M. 
Kelleher, et al., [15] and S. Winoto, et al., [16], among 
others, have been able to verify these numerical predictions 
by using hot wire anemometry, laser Doppler systems and flow 
visualization techniques to describe the motion of the secon- 
dary flow. A non-linear analysis of Gortler vortices was 
conducted by Y. Aihara [17] in 1976. 

Growing interest in the various effects of the secondary 
flow associated with Tay lor-Gort ler vortices led to more 
specific studies concerning the vortice influence on heat 
transfer in both turbulent and laminar flow. In 1955, 
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F. Kreith [4] reported that the heat transfer from a heated 
concave wall was considerably higher than that from a convex 
wall of identical curvature under similar conditions of tur- 
bulent flow. 

In the special cases of very high or very low Prandtl 
number, L. Persen [18] related the increased heat transfer 
from curved walls to the presence of Taylor-Gortler vortices . 
He determined that their overall effect was to enhance heat 
transfer through the boundary layer. 

In 1970, McCormack, et al., [19] published the results 
of their initial experimental work regarding the effect of 
secondary flow on heat transfer. From the data obtained they 
concluded that a complete theoretical explanation of the heat 
transfer effects required that all non-linear terms be re- 
tained in the flow equations. R. Kahawita and R. Meroney [20] 
considered this hypothesis and further reported that both 
the high order terms and the normal velocity components of 
the main flow were essential to any calculations at small 
wave numbers . 

In 1966, Y. Mori and Y. Uchida [5] considered the fully 
developed forced convection heat transfer between horizontal 
plates. In the results of their investigation they reported 
that as the temperature difference between the plates was 
increased a critical value developed beyond which longitudinal 
vortex rolls, similar to the Taylor-Gortler vortices, would 
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form with their axes in the streamwise direction. M. Akiyama, 
et al., [21] confirmed the presence of these longitudinal 
vortices at a critical Rayleigh number of 1708 in cases where 
the plate is heated from below. 

Although heat transfer in curved duct flow has , and will 
continue to be investigated experimentally, considerable 
emphasis has been directed recently towards heat transfer 
predictions using numerical approximation for flows in curved 
duct geometries. In 1970, K. Cheng and M. Akiyama [22] 
developed a numerical solution for forced convection heat 
transfer in laminar curved rectangular channels by using a 
point successive over-relaxation method. For similar flow 
and boundary conditions, G. Yee and J. A. C. Humphrey [23] 
used the fully elliptic forms of the transport equations to 
compose a numerical solution. However, both of these numer- 
ical solutions were limited to small aspect ratios of essen- 
tially square cross section. 

In contrast to the large volume of published numerical 
approximations to this heat transfer problem, the number of 
experimental and analytical investigations have been very 
limited. Results for a fully developed, constant wall heat 
flux, curved, square cross section channel were reported by 
Y. Mori, et al . [24], and results for larger aspect ratios 
were reported by M. Durao . [25] 
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II. INTENT OF THE STUDY 



In view of the possible engineering design implications 
of the Taylor-Gortler vortex phenomenon, the primary purpose 
of this study was to investigate the effects of the secondary 
flow on the heat transfer rates in a curved channel of rec- 
tangular cross section. 

It was expected that the presence of the Taylor-Gortler 
vortices would result in the secondary flow velocity com- 
ponents transporting the heated fluid particle from the cur- 
ved wall to the center of the flow channel. This fluid 
motion would displace the cooler fluid particle in the center 
of the channel causing it to move toward the heated curved 
wall. It was anticipated that this cyclic motion would 
measurably improve the heat transfer process for all fluid 
velocities in the laminar flow regime. Secondary investiga- 
tions were conducted to determine the additional influence 
of buoyant effects on the development of heat transfer in 
the curved channel. It was expected that the heat transfer 
rates would be greater when the buoyant forces were in the 
direction of the main flow. 

The investigation was conducted using a rectangular cross 
section channel that contained both a straight and a curved 
test section. The results obtained in the straight 
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section served as a baseline and as such were compared to 
the data obtained in the curved section. The complete 
process was accomplished for two distinct channel orienta- 
tions: Cl) buoyant forces aligned with, and (2) buoyant 

forces opposed to, the main direction of the fluid flow. 

Straight section data were compared with analytic solu- 
tions provided by P. McCuen, et al . , [26], P. M. Wors^e- 
Schmidt [27] and R. K. Shah and A. L. London. [28] The 
buoyant effects were qualitatively compared with the data 
contained in a numerical computation by R. Chilukuri and 
J. A. C. Humphrey. [29] 
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III. EXPERIMENTAL WORK 



A. DESCRIPTION OF THE APPARATUS 

A rectangular cross section channel built as described 
in reference 25 was used to achieve the objectives of this 
study. The channel was manufactured from two 0.635 cm thick 
sheets of Plexiglas separated by 0.635 cm thick spacers whose 
inside lengths served as the sides of the channel. 

As shown in Figure 2, the channel was composed of a 
straight portion 122.0 cm long, an intermediate curved 
section which formed a 180 degree arc and a final short 
straight section. The radius of curvature of the interior 
concave wall of the curved portion was 30.5 cm. The channel 
was 0.635 cm high and 25.4 cm wide with an aspect ratio of 
40 and a cross sectional area of 16.13 square centimeters. 

A cross sectional view of the channel is shown in Figure 3 . 

The working fluid was air at room temperature which 
entered the flow channel through an entrance bell constructed 
of Plexiglas . The entrance bell was designed and manufactured 
in accordance with ASME nozzle standards with an elliptic 
curve based on a ten inch major and one inch minor axis. A 
permeable cloth covered the entrance nozzle to reduce the 
introduction of dust and other foreign matter into the flow 
channel. The flow at the end of the channel exited through 
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Straight test section 
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figure 2. Schematic Illustration of Apparatus 



PLEXIGLAS mm ALUMINUM 
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f igure 3. ^ross Section of Channel at Both 'Test Sections 





an aluminum exhaust bell connected by flexible tubing to a 
Fisher and Porter Company variable area flow meter, model 
number 10A3565A. The rotameter had a 100% full scale rating 
of 0.314 cubic meters of air per minute (11.1 standard cubic 
feet of air per minute) at one atmosphere and 70° F. 

The flow of air was drawn through the channel, the flexi- 
ble tubing and the rotameter by a single electrically driven, 
Cadillac centrifugal blower, model G-12. The blower speed 
was controlled by varying the motor voltage with a General 
Radio Company Variac Autotransformer type W10MT3 . 

For the purpose of obtaining experimental heat transfer 
data, two test sections were constructed; (1) a straight 
test section, 29.2 cm in length, located in the straight 
portion of the channel downstream of the hydrodynamic entrance 
region, and (2) a curved test section 28.3 cm in length, sub- 
tending a circular arc of 53.1° in the lower quadrant of the 
curved portion of the channel. The area of the straight 
test section was 741.9 square centimeters and that of the 
curved test section was 717.7 square centimeters. 

In each test section the outer wall was detached and 

modified as shown in Figures 4 and 5. Temsheet, which is a 

carbon impregnated porous paper with the property of uniform 

/ 

electrical resistivity, was glued to the entire interior 
surface of the outer wall in each section. The flow of air 
could then be heated through the Joulean heating of the Tem- 
sheet. Since the Temsheet had an electrical resistance that 
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Figure '4. Internal Detail of the Straight Test Section 




Figure 5. Interna! Detail of the Curved Test Section 
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Thermocouples 




Insulation 
Plexiglas 
Temsheet 



283 mm 



Thermocouple 
suspended in 

the How 




was slightly temperature dependent, a precision resistor with 
an electrical resistance of 2.0236 ohms was connected in 
series with the Temsheet to facilitate calculating the 
instantaneous power supplied. A LAMBDA Regulated D.C. Power 
Supply model LK345A FM supplied electrical current to heat 
the Temsheet to desired levels. 

The variables to be measured in the investigation were: 
the flow inlet temperature at the entrance to each test 
section (T^) , the flow outlet temperature at the exit of 
each test section (Tq^) , the wall temperature of the heated 
plate in each test section (T wo ) , the wall temperature oppo- 
site the heated plate in each test section (T . ) , the flow 



exit temperature at the end of the channel (T r ^j^) , the temper- 
atures between the two insulation layers for each test section 



(T 



INS 



), and the voltages across the precision resistor (Vpp) 



and the heater (V..) . Copoer-Constantan thermocouples were 

n 

used to measure all the above temperatures. 



Liquid Crystals 

Cholesteric Liquid Crystals were employed to provide 
qualitative measurements of the surface temperature distri- 
bution on the Temsheet . The liquid crystals exhibit drama- 
tic changes in color for a small temperature differential. 
The liquid crystals would allow a reasonable qualitative 
comparison to be made between the heat transfer processes on 
each test section by observation of the color pattern 
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displayed for various temperature gradients and flow velo- 
cities. Specifically, the occurrence of multi-color stripes 
aligned in the flow direction on the curved wall would indi- 
cate regions of high and low temperatures and confirm the 
presence of Taylor-Gortler vortices. 

To select the appropriate liquid crystal system for 
application to the surface of the Temsheet, consistent 
values of flow inlet temperature, Reynolds number, Prandtl 
number and Nusselt number were assumed and an approximate 
energy balance was developed: 



Q = m C . (T_ Trrp 
x pair OUT 



- t in> 



h a pl at 



or 



AT 



T - T 
1 OUT IN 



Re Pr 
Nu 



(A ) 
c 



(A 



PL) 



where 



AT = T 



wo 



T + T 
IN uUT 



A desired ten degree Celsius difference between the 
inlet and outlet temperatures required an average wall 
temperature of forty-five degrees Celsius. Taking into 
consideration the possible losses in the system, the required 
power to achieve that average wall temperature was estimated 
to be a maximum of 40 Watts. 
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Upon completion of a series of preliminary experiments it 
was established that a combination of liquid crystals would 
be required to adequately delineate the average wall tempera- 
ture and the vortex phenomenon. The test sections were 
therefore sprayed with a mixture of three NCR Microencapsu- 
lated Liquid Crystals (R-41, R-45 and R-49) prepared in 
accordance with reference 30. Using an artist's air brush, 
approximately twenty coats of the liquid crystal mixture were 
applied directly to the Temsheet in order to ensure a sharp 
color contrast and definition within the response curve of 
each liquid crystal. It was hoped that the resultant color 
patterns could be recorded photographically. 

A total of sixty-four copper-constantan , glass insulated, 
30 gauge thermocouples were constructed and placed in specific 
locations to measure the required temperatures. All ther- 
mocouples were permanently wired to a Hewlett-Packard Data 
Acquisition System, Model 2010C, which provided an automatic 
printed record of the instantaneous data for each thermo- 
couple position. The full set of thermocouples were connected 
in series to a Kaye Instruments Ice Point Reference, Model 
K-140-4. In each test section five thermocouples were inser- 
ted between the first and second layers of the insulation 
with an additional five inserted between the second and third 
layers. Each set of five thermocouples were connected in 
parallel to read an averaged temperature value between 
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the respective insulation layers. A set of four thermo- 
couples, also connected in parallel, were inserted into the 
flow channel at each of the following locations : the en- 

trance to the channel, the exit of the straight test section, 
and the exit of the curved test section. At each location 
the thermocouples were distributed in a spanwise interval of 
5.0 cm and at a depth of 0.3 cm to allow the thermocouples 
to read a bulk temperature for the air flow. In each test 
section eight thermocouples were placed in direct contact 
with the Temsheet through very small diameter holes drilled 
in the Plexiglas above the Temsheet and the thermocouple 
beads were electrically insulated from the Temsheet with 
ENMAR Heat Resisting Glyceryl Phthalate. An additional 
eight thermocouple beads were embedded in the Plexiglas 
channel opposite those positioned in the Temsheet. Figure 
6 provides a detailed sketch of this thermocouple arrange- 
ment. The thermocouples in contact with the Temsheet pro- 
vided an average outer wall temperature (T ) and confirmed 
temperature distributions indicated by the liquid crystal 
patterns . Those in the Plexiglas channel wall provided 
an average inner wall temperature (T . ) and assisted in 
estimating the losses in the system due to radiation. All 
the thermocouples were calibrated using a ROSEMOUNT Communta- 
ting Bridge, model 920A, and a ROSEMOUNT Constant Temperature 
Bath, Model 913A. 



31 




32 



Figure 6. Location of Thermocouples in Test Sections 



The thermal insulation consisted of Armstrong ARMAFLEX 
22 Sheet Insulation which is a flexible foam plastic material. 
Three layers of insulation were used, each one with a 
thickness of 0.635 cm and covering an area Slightly larger 
than each test section. The insulation was fixed in place 
with the use of an Adhesive Heat Resistant Ventilation Duct 
Tape as shown in Figure 7. The entire length of the channel 
was sealed with GENERAL ELECTRIC Silicone Rubber sealant 
caulk . 

Two aluminum electrodes, 0.318 cm thick, were inserted 
between the Temsheet and the Plexiglas as shown in Figure 3. 
Each plate consisting of the Temsheet, the liquid crystals, 
the thermocouples and the aluminum electrodes was affixed in 
the channel using Teflon screws to provide electrical insula- 
tion . 

B. EXPERIMENTAL PROCEDURES 

Experiments were conducted for each test section following 
the same set of procedures. The data obtained in the straight 
test section provided baseline information for later curved 
section runs, either aligned with or opposed to the buoyancy 
forces. The volumetric flow rate was varied incrementally 
from 0.063 to 0.220 cubic meters of air per minute to corres- 
pond to values of Reynolds numbers in the laminar flow regime. 
Data was recorded and calculations were performed for each 
flow rate increment. 
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Figure 7. General Arrangement of Apparatus and Instrumentation 





The appropriate electrical power required to heat the 
Temsheet to a temperature level that would initiate the 
liquid crystal response was determined through preliminary 
experiments . 

Estimates of the time required for the experiment to 
reach steady state, decisions on what data would be recorded, 
and frequencies at which that data would be obtained were also 
based on preliminary runs. 

The determination as to what constituted a steady-state 
condition required careful consideration in light of the in- 
fluence that the multiple input variables could have on the 
results. The final criteria for steady-state was based on 
three input variables; the flow rate percentage, the outside 
wall temperature (Temsheet), and the heater voltage. When 
these input conditions had a variation of less than two 
percent over a ten minute interval, the experiment was consi- 
dered to have reached steady-state. Under these guidelines 
the assumptions of steady flow and constant heat flux from 
the upper wall were valid. 

The time required to reach steady- state was found to be 
between two and two and one-half hours. Accordingly, all 
runs were two and one-half hours in duration. Data was 
recorded at one hour intervals for the first two hours and at 
ten minute intervals thereafter. 
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The instantaneous values for all the temperatures were 
automatically printed by the digital recorder in milivolts . 

A standard thermocouple table was used in conjunction with 
previously obtained thermocouple calibration data to convert 
the milivolts to degrees Celsius. The volumetric flow rate 
was read directly from the rotameter. The precision resis- 
tor voltage (V pR ) and the heater voltage (V^) were recorded 
in volts and the instantaneous voltage supplied (0 ) was 
calculated using the following relationship: 

CL 



V V 
PR H 



R 



PR 



where R pR was 2.0236 ohms, the measured electrical resistance 
of the precision resistor in the circuit. 

The steady state response of the liquid crystals was 
recorded photographically using a NIKOMAT (FTN) 35mm single 
reflex camera with a 85 - 205 mm zoom lens. The film was 
25 ASA Kodachrome with the shutter speed and lighting con- 
trolled manually. 

All data reduction was accomplished by desk-top computer 
programs and stored routines. 
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IV. PRESENTATION OF DATA 



A. ANALYSIS 

The flow examined in this study was laminar and steady, 
thus in a general sense, plane Poiseuille flow. The Temsheet 
which possessed the property of uniform electrical resistivity 
approximated a constant axial and peripheral wall heat flux. 
Temperature measurements on the Plexiglas' wall opposite the 
heated surface indicated negligible heat loss through that 
wall, i.e., adiabatic. With an aspect ratio of 40 the chan- 
nel was approximated by infinite parallel plates. The 
straight portion of the flow channel, upstream from the 
straight test section, was of sufficient length to ensure 
that the flow was hydrodynamically developed for all flow 
velocities in the laminar range. For the purposes of this 
study, the flow was considered laminar when the value of the 
Reynolds number was less than 2000. Based on the above 
assumptions, the experimental configuration was then modeled 
as laminar forced convection between parallel plates with 
hydrodynamically developed and thermally developing flow. 

The boundary conditions are fundamental boundary conditions 
of the second kind with one wall at a specified axial and 
peripheral heat flux and the other wall insulated. [28] 
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The appropriate variables, corresponding equations and solu- 
tion procedures are contained in the following paragraphs. 

The average heat transfer coefficient from the heated 
wall to the flow of air was defined by a form of Newton's 
law of cooling: 



Q . = h A DT AT 

x air PL 

In this expression, Q • was the heat convected to the air, 
h was the average heat transfer coefficient, Ap^ was the 
surface area of the heated Temsheet, and AT was the tempera- 
ture difference between the average wall temperature (T ) 
and the bulk temperature of the fluid (Tg). 

The actual heat convected to the air (Q . ) was then 

cLUT* 

independently calculated with the equation: 



Q 



air 



= mC 



(T 



-T t „) 



pair OUT IN 



where C . was the specific heat of air at constant pressure 
pair r r 

and m was the mass flow rate of the air. The mass flow rate 
was the product of the volumetric flow rate (q) and the local 
flow density in the rotameter (p). The volumetric flow rate 
was corrected to standard air conditions for the rotameter 
and the local flow density was calculated assuming perfect 
gas behavior at constant pressure and using the flow exit 
temperature • 
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The average Nusselt number was determined from the aver- 
age heat transfer coefficient using the following expression: 



Nu 



hd 

K • 



air 

where d represented the height of the channel and K . was 
the thermal conductivity of the air evaluated at the bulk 
temperature (T n ) . 

The Reynolds number was calculated for all the experimen- 
tal runs and was defined as: 



R pud = md 
r laA 0 

where d was the height of the channel, y was the dynamic 
viscosity of the air at the flow exit temperature and A 
was the cross sectional area of the channel. 

For all experiments in the curved test section, the Dean 
number was calculated as : 



De 



Re 




In this equation d was again the height of the channel and 
Ri was the radius of curvature of the inner convex wall. 

To verify that Q a j_ r J which was indirectly measured by the 
temperature rise across the channel, was actually the con- 
vection heat transfer from the heated Temsheet surface , an 
estimate of the radiation heat transfer was made for each run. 
The average temperature of the lower plate was measured by 
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eight thermocouples and the radiated heat transfer (Q- ) was 
calculated with the Stefan-Boltzman law of radiation: 



<*r * 



4_ t - 4 ) 

WO Wl 



where a is the Stefan-Boltzman constant and represents 
the total resistance to radiation heat transfer between the 
plates as given by: 



R R 





wo wi 



1 ] 



In this form e represented the emissivity of the heated 
wo 

Temsheet surface and e • was the emissivity of the adiabatic 

Wl 

Plexiglas surface . 

Appendix B contains a representative heat balance and 
sample calculations for one of the experimental runs. 



B. RESULTS 

Utilizing the expressions described earlier, the experi- 
mental data for buoyant forces aligned with the main direction 
of the fluid flow was reduced in accordance with the sample 
calculation procedure shown in Appendix B. Results for the 
major variables involved are provided in tabular and graphical 
form on the following pages. Table I and II show computed 
values for the straight and the curved test sections respec- 
tively. Figure 8 contains a plot comparing the average 
Nusselt number versus the Reynolds number for both test 
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sections. Representative uncertainty bands have been indi- 
cated for all plotted experimental data and a complete error 
analysis is provided in Appendix A. 
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TABLE I 



SUMMARY OF RESULTS FOR STRAIGHT TEST SECTION 



Re 


h(W/m 2 °C) 


Nu 


989 


22.51 


5 .321 


876 


22 . 21 


5 . 298 


826 


21 . 79 


5 .198 


751 


21.75 


5 .197 


689 


20.72 


4.941 


632 


20.95 


4.897 


560 


19 . 20 


4.559 


488 


18 . 79 


4.482 


438 


17 .75 


4.234 


363 


16.99 


4.053 


313 


16.24 


3 .874 


263 


14 . 30 


3 . 512 
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TABLE II 



SUMMARY OF RESULTS FOR CURVED TEST SECTION 



Re 


Dc 


h( W/ni °C ) 


Nu 


990 


144 


27.81 


6 .694 


865 


126 


27 . 74 


6 .617 


803 


117 


25.06 


5.977 


741 


108 


23.42 


5.588 


678 


99 


22 .57 


5 .384 


654 


95 


22.43 


5 .351 


616 


90 


22.07 


5.265 


579 


84 


21.67 


5 .169 


554 


81 


21.34 


5 .091 


523 


76 


20.58 


4.875 


492 


72 


20 . 28 


4.838 


429 


63 


20 . 20 


4.819 


404 


59 


20 .31 


4.844 


367 


54 


18.95 


4 . 515 


342 


50 


17.76 


4 .236 


306 


45 


15.76 


3 . 759 


274 


40 


14.50 


3 . 460 


244 


36 


13.12 


3.130 
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In all cases the results indicated an increase in the 
heat transfer rates as the flow velocity, i.e., Reynolds 
number, was increased. However, when comparing the average 
Nusselt numbers of the straight and curved sections two 
distinct regions became evident. At the highest values of 
Reynolds number, which corresponded to the highest volu- 
metric flow rates, it was observed that the values of the 
average Nusselt number for the curved test section were 
greater than those in the straight test section. At the low 
and intermediate values of Reynolds number, the data was 
less conclusive with regards to the comparative values of 
the average Nusselt number in each test section. This region, 
between Reynolds numbers of 200 and 400, corresponded to 
values of the Dean number between 30 and 60. As Dean discussed 
in his work concerning flow instability in a curved channel 
[12], a small disturbance would not be amplified below the 
critical Dean value of 36. Above this critical value, his 
linearized prediction indicated the disturbance would increas- 
ingly amplify as flow velocity increased. 

Consequently, it was concluded that the Taylor-Gortler 
vortices are only initially developing in the low to inter- 
mediate Reynolds number region and have not yet reached suf- 
ficient proportion to have any significant effect on the 
heat transfer rates in the curved section. 
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The actual presence and onset of the Taylor-Gortler 
vortices was independently verified through the observation 
of the liquid crystal isotherm distribution on the Temsheet 
surface in the curved test section. At the smallest values 
of Reynolds number, the liquid crystals were in essentially 
straight bands that were perpendicular to the flow direction 
(Figure 9). A similar pattern was observed in the straight 
test section for the entire experimental range of Reynolds 
numbers. This pattern was thought to indicate that the flow 
was in essence two-dimensional. As the value of the Reynolds 
number was increased, liquid crystal bands orientated in the 
streamwise directions became visible at the trailing edge of 
the curved test section. At the onset of these streamwise 
bands the liquid crystal bands at the leading edge remained 
orientated perpendicular to the flow direction. As the 
Reynolds number was further increased, the bands which were 
orientated in the streamwise direction eventually became 
evident over the entire surface of the curved test section 
(Figure 10). These streamwise bands would indicate a transi- 
tion from the initial, stable, two-dimensional flow to a less 
stable, three-dimensional flow triggered by the continuing 
de velopment of the Taylor-Gortler vortices . 

The simultaneous occucrencs of the streamwise bands, 
explained by the presence of the longitudinal Taylor-Gortler 
vortices, and the increase in the average Nusselt number 
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Figure 9. Distribution of Liquid Crystal 

Color Bands in Curved Test Section 
for Re = 2 5 0 (De = 35) 
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Figure 10 . Distribution of Liquid Crystal 

Color Bands in Curved Test Section 
for Re=800 (De=117) 
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justified the conclusion that the vortices did in fact en- 
hance heat transfer from the curved test section. 

Comparison of the average Nusselt numbers for the 
straight and curved test section indicated that values for 

the ratio of Nu to Nu could be separated into two cate- 

c s 

gories. For Reynolds numbers between 400 and 800, the heat 

transfer rate was increased by an average of eleven percent. 

For Reynolds numbers above 800, the heat transfer rate was 

increased by more than twenty-five percent. In an earlier 

investigation of the boundary layer along a concave wall, 

Kreith [4], reported that heat transfer rates along concave 

walls increased from twenty-five percent to sixty percent 

for Reynolds numbers, based on the hydraulic diameter, 

4 6 

between 10 and 10 . This would suggest that the two data 
points for Reynolds numbers greater than 800 in the curved 
test section could represent flow that has tripped and 
become turbulent, or at least entered a laminar-turbulent 
transition. One possible explanation could be that the 
Taylor-Gortler vortices have started to break down at this 

point and have therefore initiated transition to the turbulent 

« 

flow regime. 

Experimental data obtained for the curved test section 
with the buoyant forces opposed to the main direction of the 
flow indicated a minimal difference in heat transfer rates 
from the curved channel. This small decrease was only 
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evident at very low volumetric flow rates and became negli- 
gible for Reynolds numbers greater than 500. 
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V. DISCUSSIONS AND CONCLUSIONS 



Investigation of the measured wall and air flow tempera- 
tures indicated there was negligible difference between the 
unheated wall temperature and the fluid bulk temperature. 

In addition, energy balances such as the one detailed in 
Appendix B provided evidence that the radiation heat transfer 
from the unheated wall was minimal. Based on this information, 
it was assumed that the heat convected to the air flow was 
solely from the heated wall. 

As previously discussed, the high aspect ratio of the 
channel provided the basis from which the initial assumption 
was made to model the experimental apparatus configuration as 
infinite parallel plates. The experimental data was found to 
substantiate this assumption as the computed values of the 
average Nusselt number in the straight section asymptotically 
approached the theoretical value of 5.385 for parallel plates 
with one wall at a constant heat flux and the opposite wall 
adiabatic. [28] 

With the formulation of the problem complete, the experi- 
mental data could now be compared with analytical solutions 
of the same problem; i.e., laminar flow between infinite 
parallel plates with one wall at constant heat flux and the 
opposite wall adiabatic. In particular, the experimental 
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results were compared with two analytic studies; one con- 
ducted by McCuen, et al. [26], for infinite parallel plates, 
and the other conducted by Wors^e-Schmidt [27] for tubes 
and annular passages. In the case of the annular passages 
the results were compared to the numerical predictions for 
an annulus with an infinite inner and outer radius of curva- 
ture, thereby approximating parallel plates. 

Before comparing the experimental data with the studies 
mentioned above, it would be advantageous to address the 
initial assumptions and review the nomenclature common to 
both those studies. First, the following restrictions are 
applied in the mathematical formulation of the problem: 

1. the velocity profile is fully established 

2. the fluid transport properties and density are 
assumed constant 

3. axial heat conduction is negligible 

4. viscous energy dissipation is negligible 

5. conditions are invariant with time (steady state) 

All of these restrictions have been accounted for in the 
design of the experimental configuration and/or verified in 
the experimental process. Next, a dimensionless axial length 
coordinate was defined as follows: 

x* = x 

Oh Pr Re 
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where x was the length of the straight test section (29.2 cm), 
Dh was the calculated hydraulic diameter (1.239 cm), and Pr 
was the Prandtl number (0.70). In these studies, the 
Reynolds number (Re) was calculated on the basis of the 
hydraulic diameter. The values of x* indicated that the 
development of the experiments would most likely occur in 
the thermal entrance region. Analysis of the temperature 
distribution on the heated wall versus the axial length 
verified the flow was indeed thermally developing. 

In comparing the experimental results to those reported 
by McCuen, et al . , it was necessary to integrate their 
expression of the local Nusselt number for small values of 
x* ; 

Nu 5 777 — — 

0.67095 x* 1/ 3 -2( 1 + 7=r^)x* 

^wo 

in order to obtain the following expression for the average 
Nusselt number. 



Nu = 



--3 



Q w , 

4(l+fpi) x* 

^WO 



-i _ r 2 /■ , j_ ^wi ^ ... 2/ 3 

ln “ 0.67095 (1 Q )x “ - 1 

wo 



In this expression Q . would represent the heat convected to 
the air from the inner Plexiglas wall and Q wq would be the 
heat convected from the outer Temsheet wall. The average 
Nusselt number was then calculated in accordance with the 
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experimental model where Q . was equal to zero and Q ^ was 
r W1 wo 

synonymous with Q a £ r * 

The experimental data obtained for the straight section 
was then recalculated as a function of the hydraulic diameter 
and both the experimental and analytical results were plotted 
in Figure 11. Note that the analytic expression for the 
average Nusselt number was only plotted for values of the 
Reynolds number that corresponded to x* less than 0.042 as 
this represented the maximum value for which the analytic 
solutions would be valid. Neglecting those experimental data 
points that plotted outside this valid range, there was good 
agreement between the experimental and the analytic results. 

A similar process was followed to compare the experimental 
results with those reported by Wars# e-Schmidt . Again the 
analytic values for the local Nusselt number were integrated 
to obtain analytic values for the average Nusselt number. 

In this case the experimental and analytic values of average 
Nusselt number were plotted versus the non-dimensional axial 
length coordinate x* . The results are shown in Figure 12. 
Similarly, there was good agreement between the experimental 
and the analytic results . 

It is noted that in both cases where the analytic and 
experimental results were compared, the experimental data 
plotted above the analytic predictions. This difference can 
be explained, at least in part, by the limitations that are 
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Figure 11. Comparison of Present Results with the Study by McCuen, 







II 

* 

X 



A 

Q 



D 



12 



56 



Figure 12. Comparison of Present Results with the Study by Wors^e-Schm i dt 



inherent in any experimental configuration. In this case, 
the presence of side wall effects and the inability to 
totally account for all the heat transfer processes and/or 
losses combine to produce a slightly inflated value of the 
heat transfer rate. Consequently, the average Nusselt number 
is also slightly inflated. 

Irrespective of the slight differences locally, the over- 
all trends were conclusive. The heat transfer rates were 
shown to increase with increasing values of Reynold number. 

The ratio of N to N indicated that there was a measureable 
uc us 

increase in the heat transfer rate from the curved section 
once the Taylor-Gortler vortices had developed to sufficient 
strength. Although much less conclusive from this experi- 
mental data, it appeared the the Taylor-Gortler vortices 
were additionally important in the transition from laminar 
to turbulent flow. 

The experimental data from the curved test section when 
the buoyant forces were opposed to the main flow direction 
was essentially inconclusive. While the trend was qualita- 
tively what had been expected, that is, the average Nusselt 
number for the buoyant forces aligned with the flow was 
greater than that for the buoyant forces opposed to the flow; 
the difference was insignificant in light of the experimental 
uncertainties . 
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VI. RECOMMENDATIONS 



While the amount of work remaining in the experimental 
study of heat transfer in a curved channel has been reduced 
for laminar flows, a large amount still remains to be 
accomplished. Additional experiments should be conducted 
to refine the correlation for the average Nusselt number as 
a function of the Dean number. Further work should certainly 
be considered in the region where the laminar flow begins its 
transition to turbulent flow. Detailed study in this area 
could provide important data to determine the manner of 
relationship that exists between Taylor-Gortler vortices and 
laminar-turbulent transition. With this work it will be 
necessary to design and construct a larger experimental 
configuration capable of higher volumetric flow rates and 
increased AT through the test sections. 

Finally it would be advantageous to do more experimental 
work in the general area of liquid crystal thermography, and 
more specifically, the photography of these crystals and their 
distribution . 
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APPENDIX A 



ERROR ANALYSIS 



The uncertainties for the major variables in the experi- 
ments were calculated in accordance with the method described 
by S. Kline and F. McClintoch [31]. The estimates of the 
uncertainty in the measured quantities were made conserva- 
tively so that there was considerable confidence in the cal- 
culated uncertainties. As an example, the calculation of the 
uncertainty for the Reynolds number is given below. The 
Reynolds number was defined by the equation: 



and the uncertainty was calculated as: 



The uncertainty in the mass flow rate was determined as 
follows : 



The uncertainty in the pressure, the gas constant, the tem- 
perature, the reading of the volumetric flow rate, the 
height of the channel, the dynamic viscosity and the cross 
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sectional area of the channel were obtained from estimates as 
follows : 



dP 

P 



dd 
d " 



0 . 00048 , 



0 . 02003 , 



dR 

R 



dy 



air 



air 



0.00002 , =£- 
= 0.00223, 



= 0.00056, 



and 



d A 

c _ 

A 



c 



^ = 0.020 

0.02054. 



The uncertainty in the mass flow rate was: 



y (0. 00 048 ) 2 +(0. 00002) 2 +(0. 00056) 2 +(0. 02000) 2 ': 

and the uncertainty in the Reynolds number could then be 
calculated as: 




dRe 

Re 



■V 



(0.02001) +(0.020 03 )+(0.00223)+(0.02054)=0 



The values of the uncertainties for all variables ar 



listed below: 



Quantity 



A 

c 

A PL 

^t>air 

De 



K . 



air 

Nu 



Uncertainty 

0.02054 
0.00266 
0 .00415 
0.04042 
0 .02149 
0.00038 
0 .02938 



0.02001 



03505 
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Quantity 



Uncertainty 



m 


0 . 02001 


Q 


0 . 02000 


o . 


0.03548 


N air 


R 


0 . 03505 


e 


T 

3 


0.00482 


T 

IN 


0 . 00257 


T 


0 .00950 


OUT 


T 


0 .00583 


w 


T T 


0 .02900 


OUT- IN 


at =V t b 


0 . 00736 


y 


0.00223 


p 


0.00074 
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APPENDIX B 



SAMPLE CALCULATIONS 

Figure 13 provides a representation of the control volume 
used in the energy balance on the straight test section. The 
sample calculations that follow indicate the procedure used 
to obtain the heat transfer components required for the 
energy balance and eventually the values of the average 
Nusselt number as a function of the Reynolds number for the 
air flow. A similar energy balance and procedure was used 
for the curved test section. 

ENERGY BALANCE 




Figure 13. Energy Balance in the Straight Test Section 



SAMPLE CALCULATIONS 



A . DATA 



t room 


= 


23.0° 


C 


Q 


= 


45% 




V PR 


= 


2 .583 


V 


V 

U 


= 


30.008 V 


r pr 


= 


2.024 


ft 


t in 


= 


0 .892 


mV 


t insi 


= 


1.597 


mV 


T INS2 


= 


1.339 


mV 


Twol 


= 


1.853 


mV 


Two 2 


= 


1.843 


mV 


Two 3 


r 


2 . 010 


mV 


Two4 


= 


2.014 


mV 


Two 5 


= 


2.144 


mV 


Two6 


= 


2.141 


mV 


Two7 


- 


2.274 


mV 


Two 8 


= 


2.239 


mV 


Twil 


= 


1.079 


mV 


Twi 2 


= 


1.083 


mV 


Twi3 


= 


1.130 


mV 


Twi4 


= 


1.131 


mV 


Twi 5 


= 


1.183 


mV 


Twi 6 


- 


1.172 


mV 



= 23.07 °C 

= 39.54 °C 

= 33.45 °C 



Two = 50.47 °C 



Twi = 28.62 °C 
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Twi7 


= 


1.219 mV 




Twi8 


= 


1.216 mV 




t out 


r 


1.383 mV = 


34.09 °C 


t exit 


= 


0.934 mV = 


23.53 °C 


a pl 


= 


0.0742 m 2 




A 

c 


= 


0.0016 m 2 




d 


= 


0.00635 m 




K air 




0.02662 w/m 


°C 


Cp air 


= 


1.0057 KJ/Kg °C 


p air 


= 


1.983 x 10" 5 


Kg/m* sec 


k ins 




4.18 x 10" 2 


W/m °C 


ax ins 


= 


0.00635 m 




£ 

WO 


= 


0.70 




£ . 
W1 




0.40 




F 

WO-Wl 


= 


1.0 




a 


= 


5.669 x 10" 8 


W/m 2 °K 



B. 



TEMPERATURE CALCULATIONS 

1 . Temperature Difference in Insulation 



2 . 



3 . 



AT INS = T INS1~ T INS2 39.54-33.45 
Bulk Temperature (Tg) 

~ _ T IN + T 0UT _ 23.07 + 34.09 _ 

B " 2 2 

Mean Temperature Difference ( AT ) 

AT = T - T = 50.47 -28.58 = 21.89 
wo B 



(AT INS ) 
6.09 °C 



28,58 °C 



°C 



64 



C. POWER CALCULATIONS 



1. Power Supplied (Qp) 

„ _ V PR V H _ (2.583) (30.008) 

^P " R pR (2. 024) 



38 . 30 W 



2. Heat Lost Through Outer Plate (Q^) 

AT INS _ 6.09 

^Lo ' AX TM0 /K_. T0 .A t5T " (0.00635)/ ((4. 18xl0- 2 ) ( 0.0742)) 
INS INS PL 



2 . 97 W 



3. Heat Radiated (Q ) 



a. Radiation Resistance (R R ) 



1-e 



R T 



wo 



1-e • 

wi 



■ A PL £ wo 



A„ t F . A_ t e . 

PL wo-wi PL wi 



i c^. ♦ . i: = ijpi 

a pl 



A nt e e . 

PL wo wi 



b. Heat Radiated (Qr) 



Qr = 



4 4 

a(T -T . ) 

WO Wl 

r r 



T = °K 



( 5. 6 6 9x10 ~ 3 ) (0.0742) (323. 63^-301. 78 4 ) 

2.929 



= 3.843 W 
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4, Heat Convected to Air (Q a ^ r ) 



a. Corrected Volumetric Flow Rate (q) 

/ 



q = Q/0 = Q 



T (°R) 

EXIT V ; 

T f O p A 

1 STANDARD ^ 



( 0 . 314) (0.45) / 534.04 

60 / 529.69 

3 

= 0.0024 m / sec 

b. Density (p) 

? 1.013xl0 5 

? ~ R T exi “ " (287 M 23. 53 + 273.16)' 



= 1.189 Kg/m 3 



c. Mass Flow Rate (m) 

m = qp = (0.0024) (1.189) = 0.0028 Kg/sec 

d. Heat Convected to Air (Q a ^ r ) 

^air = m Cpair (T 0UT" T IN ) 

= (0.0028)(1.0057)(34. 09-23. 07) 

= 0.03103 KJ/sec 

= 31.03 W 
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D. AVERAGE HEAT TRANSFER COEFFICIENT (h) 



h = 



^air 
ApL AT 



31.03 

(0 .0742) ( 21. 39 ) 



= 19.11 W/H 2 °C 



E. AVERAGE NUSSELT NUMBER (Nu) 



Nu = 



hd 

Kalr 



(19 .11) (0 .00635) 
(0.02662) 



4.559 



F. REYNOLDS NUMBER (Re) 

_ m d ( 0.0028)(0. 00635) 

Re = r— = = 

y air c (i:983xl0“°)(0.0016) 

= 560.39 
= 560 
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